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Abstract 

Layered perovskites Sr2lr04 and Ba2lr04 are regarded as the key materials for understanding 
the properties of magnetic relativistic insulators, mediated by the strong spin-orbit (SO) coupling. 
One of the most fundamental issues is to which extent these properties can be described by the 
superexchange (SE) model, formulated in the limit of the large Coulomb repulsion for some ap¬ 
propriately selected pseudospin states, and whether these materials themselves can be classified 
as Mott insulators. In the present work we address these issues by deriving the relevant models 
and extracting parameters of these models from the first-principles electronic structure calculations 
with the SO coupling. First, we construct the effective Hubbard-type model for the magnetically 
active t2g bands, by recasting the problem in the language of localized Wannier orbitals. Then, 
we map the so obtained electron model onto the pseudospin model by applying the theory of SE 
interactions, which is based on the second-order perturbation theory with respect to the transfer 
integrals. We discuss the microscopic origin of anisotropic SE interactions, inherent to the compass 
Heisenberg model, and the appearance of the antisymmetric Dzyaloshinskii-Moriya term, associ¬ 
ated with the additional rotation of the IrOg octahedra in Sr2lr04. In order to solve the pseudospin 
Hamiltonian problem and evaluate the Neel temperature (Tat), we employ the non-linear sigma 
model. We have found that, while for Sr2lr04 our value of T^r agrees with the experimental data, 
for Ba2lr04 it is overestimated by a factor two. We argue that this discrepancy is related to 
limitations of the SE model: while for more localized t2g states in Sr2lr04 it works reasonably 
well, the higher-order terms in the perturbation theory expansion play a more important role in 
the more “itinerant” Ba2lr04, giving rise to the new type of isotropic and anisotropic exchange 
interactions, which are not captured by the SE model. This conclusion is supported by unrestricted 
Hartree-Fock calculations for the same electron model, where in the case of Ba2lr04, already on 
the mean-field level, we were able to reproduce the experimentally observed magnetic ground state, 
while for Sr2lr04 the main results are essentially the same as in the SE model. 
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I. INTRODUCTION 


5d transition-metal oxides have attracted a considerable attention as a new paradigm of 
relativistic magnetic materials, whose properties are largely influenced by the strong spin- 
orbit (SO) coupling, leading to the experimental realization and a number of theoretical 
proposals for such fascinating phenomena as SO interaction assisted Mott state in Sr 2 lr 04 ,-'- 


3. y and y). 


possible existence of 


spin-liquid state in Pr 2 lr 207 (Ref. U) and Na4lr308 (Refs, 
topological semimetallic states in pyrochlore iridates,- and unusual magnetic ordering in the 
honeycomb compounds Na 2 lr 03 and Li 2 lr 03 ,-^— which may be relevant to the Kitaev model 
of bond-dependent anisotropic magnetic coupling.— 


In this respect, a lot of attention is being focused on the properties of tetravalent irid¬ 
ium oxides (or iridates), originating from the 5/6-£lled Ir t 2 g band, located near the Fermi 
level. The strong SO interaction splits the atomic t 2 g states into the fully occupied four¬ 
fold degenerate j = 3/2 states and twofold (Kramer’s) degenerate j = 1/2 states, which 
accommodate one electron. In this sense, there is a clear analogy with the spin-1/2 systems 
and the problem of interatomic exchange interactions can be formulated in terms of some 
appropriately selected pseudospin states. In solids, each group of states form the bands, 
which can, however, overlap with each other. Moreover, since j is the band quantum num¬ 
ber in solids, there is always a finite hybridization between these two groups of relativistic 
states. The j = 1/2 electrons experience the on-site Coulomb repulsion and can polarize 
the occupied j = 3/2 shell via the intraatomic exchange interactions. Moreover, the precise 
division of the t 2 g states into the j = 3/2 and j = 1/2 ones depends on the crystal-field 
splitting, which is typically smaller than the SO coupling. These are the main ingredients, 
which predetermine the low-energy electronic properties of iridates. 


The layered perovskites, Sr 2 lr 04 and Ba 2 lr 04 , are typically regarded as the key mate¬ 
rials for revealing the basic microscopic mechanisms, which can operate in iridates. They 
are also used as the benchmark materials for testing the new theoretical models. In this 
respect, the first and one of the most successful theoretical models for iridates was based 
on the theory of superexchange (SE) interactions, which is valid in the limit of large on-site 
Coulomb repulsion and treats the transfer integrals between the relativistic pseudospin states 
in the second order of perturbation theory.-ii^ This model indeed reveals a rich and very 
interesting physics, including the bond dependence of the anisotropic exchange couplings 
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and emergence of large antisymmetric Dzyaloshinskii-Moriya (DM) interactions when the 
inversion symmetry is broken by the anti-phase rotations of the IrOg octahedra in Sr 2 lr 04 . 

At the same time, there was always a qnestion abont how far one can go in applying the 
SE model for the real iridates. For the layered iridates, this point was risen in Ref. [ij, where, 
nsing the dynamical mean-held theory (DMFT), the anthors of this work have argued that 
the behavior of both Sr 2 lr 04 and Ba 2 lr 04 retain many aspects of Slater insulators, whose 
insulating properties are closely related to the existence of the long-range antiferromagnetic 
(AFM) order. The problem reemerged again recently after the experimental discovery of the 
magnetic ground state structure of Ba 2 lr 04 ,— which cannot be described by the SE model, at 
least on the mean-held level.— Thus, the questions is whether this problem should be resolved 
by considering the quantum huctuation ehects, but within the SE model,— or revising the 
SE model itself by including to it some new higher-order terms in the perturbation theory 
expansion. The answer to this question is not obvious, because in the SE formulation, 
the ehects of the SO coupling are included to the transfer integral. Therefore, the second- 
order perturbation theory with respect to the transfer integrals automatically means that 
it treat the SO coupling also only up to the second order. If the SO interaction is large 
(as in iridates), it can be rather crude approximation, because it does not take into account 
several important ehects, such as the in-plane anisotropy in the uniaxial systems, which 
may be relevant to the experimentally observed behavior of Sr 2 lr 04 and Ba 2 lr 04 . Another 
interesting point is the value of Neel temperature (Tat), which is remarkably close in both 
considered systems (about 240 K), and whether this fact can be rationalized on the basis of 
SE theory. 

The main purpose of this work is to critically reexamine abilities of the SE theory for the 
layered iridates. This is certainly not the hrst attempt to derive parameters of interatomic 
exchange interactions using the theory of SE interactions and the basic ideas of this method 
in the case of the strong SO coupling are well understood today, at least for the models.— 
Nevertheless, besides the SO coupling, the behavior of interatomic exchange interactions 
strongly depends on the number of adjustable parameters, used in the model Hamiltonians, 
such as the on-site Coulomb and exchange interactions, tetragonal crystal-held splitting, and 
the matrices of transfer integrals. Therefore, we believe that, in the process of derivation of 
the pseudospin model, it is very important to reduce the number of possible ambiguities by 
sticking as much as possible to the hrst-principles electronic structure calculations. 
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The rest of the article is organizes as follows. In Sec. [TTl we will briefly discnss the main 
differences of the crystal and electronic strnctnre of Ba 2 lr 04 and Sr 2 lr 04 . Then, in Sec. Illll we 
will explain onr method of the constrnction of the effective low-energy electron model on the 
basis of hrst-principles electronic strnctnre calcnlations with the SO conpling. This model 
will be fnrther nsed in Sec. IIVI as the staring point for the derivation of the SE Hamiltonian 
in the basis of psendospin states. In Sec. 0 we will discnss resnlts of onr calcnlations of the 
SE interactions and their implications to the magnetic properties of Ba 2 lr 04 and Sr 2 lr 04 
nsing the non-linear sigma model. Then, in Sec. lYD we will provide a detailed comparison 
with the resnlts of nnrestricted Hartree-Fock (HE) calcnlations, which do not rely on the 
pertnrbation theory, and argne that while for Sr 2 lr 04 the SE theory works reasonably well, 
for Ba 2 lr 04 it misses several important interactions, which are nonetheless captnred by the 
HE calcnlations. Finally, in Sec. I VI II we will give a brief snmmary of onr work. Details of 
derivation of the non-linear sigma model for the compass Heisenberg model will be given in 
the Appendix. 


II. MAIN DETAILS OF CRYSTAL AND ELECTRONIC STRUCTURE 


In this work we nse the experimental structure parameters, reported in Refs.l2l|andl22 (at 
13 K) for Ba 2 lr 04 and Sr 2 lr 04 , respectively. According to these data, Ba 2 lr 04 crystallizes in 
the undistorted tetragonal Am/mmm structure with the Ir-O-Ir angles in the xy plane being 
equal to 180°. Sr 2 lr 04 exhibits the additional rotation of IrOe octahedra (the space group 
Mi/acd), which leads to the deformation of the Ir-O-Ir bonds in the xy plane (see Fig. [1]). 
Depending on the Ir site, this rotation can be either clockwise (+ 0 ) or counterclockwise 
(— 0 ). The experimental value of the angle 0 is 12°.— 

The corresponding electronic structure in the local-density approximation (LDA) with the 
SO coupling is displayed in Figs. |2] and |3] for Ba 2 lr 04 and Sr 2 lr 04 , respectively. In this work 
we will focus on the behavior of magnetically active Ir t 2 g bands, located near the Fermi level 
and separated relatively well from the rest of the spectrum. There are two main differences 
between Ba 2 lr 04 and Sr 2 lr 04 : (i) The Ir t 2 g band is narrower in Sr 2 lr 04 (the bandwidth is 
about about 3 and 3.5 eV in Sr 2 lr 04 and Ba 2 lr 04 , respectively). This is generally consistent 
with the additional distortion in Sr 2 lr 04 , which leads to the deformation of the Ir-O-Ir bonds; 
(ii) The Ba 5d band in Ba 2 lr 04 , which strongly hybridizes and, therefore, has a large weight 
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FIG. 1 . (Color online) (Left) Rotations of IrOg octahedra in the xy plane of Sr2lr04. The Ir atoms 
are indicated by the big (red) spheres and the oxygen atoms are indicated by the small (yellow) 
spheres. The sites aronnd which the octahedra are rotated clockwise (+</>) and counterclockwise 
(—(/>) are denoted as 1 and 2 , respectively. (Right) The directions of axes in the lAi/acd {x,y) and 
4m/mmm {x',y') coordinate frames. 

of the Ir 5 d states, is much closer to the Fermi energy than the Sr Ad band in Sr2lr04. This 
is mainly related to the larger Ba 5 d bandwidth, due to the less distorted crystal structure 
as well as the relativistic effects.— 

In Ba2lr04, the relativistic j = 3/2 and j = 1/2 subbands are separated by the direct 
gap, which allows us to construct both 6 - and 2 -orbitals models (for the entire t2g bands 
and j = 1/2 subband, respectively). In Sr2lr04, due to the additional mixing between the 
j = 3/2 and j = 1/2 states, caused by the lAi/acd distortion, such separation does not take 
place. Therefore, for Sr2lr04, we will focus only on the 6 -orbital model. 

III. CONSTRUCTION OF EFFECTIVE LOW-ENERGY ELECTRON MODEL 

In this section we will discuss the construction of the low-energy electron model, starting 
from the LDA band structure with the SO interaction. For practical calculations, we use 
the linear muffin-tin orbital (LMTO) method in the nearly orthogonal representation.— The 
model itself has the following form: 

nei = EZ (1) 

ij ajS i ol ^')5 
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FIG. 2 . (Color online) Electronic structure of Ba2lr04 in LDA with the SO coupling, (a) Total 
and partial densities of states. The shaded area shows the contributions of the Ir 5 d states. The 
positions of the main bands are indicated by symbols, (b) Band dispersion near the Fermi level, 
as obtained for the full LDA Hamiltonian in comparison with the six- and two-orbital models. 
The high-symmetry points of the Brillouin zone are denoted as F = ( 0 , 0 , 0 ), X = (vr/a, vr/a, 0 ), 
N = (tt/o, 0 , tt/c), P = (tt/o, vr/a, vr/c), and Z = ( 0 , 0 , 27 r/c). The Fermi level is at zero energy 
(shown by dot-dashed line). 

where and are, respectively, the creation and annihilation operators of an electron 
on the Wannier orbitals Wia, centered at the Ir site i and specihed by the index a = (m, s), 
which numbers Kramer’s doublets m= 1 , 2 , or 3 (an analog of orbital indices without SO 
coupling) and the states s= 1 or 2 within each such doublet (an analog of spin indices 
without SO coupling). 

First, we construct the Wannier functions for the magnetically active bands, using the 
projector-operator technique.—^— We consider the 6 -orbital model for the both Ba2lr04 
and Sr2lr04. Moreover, for Ba2lr04 it is also possible to construct the 2 -orbital model, 
by considering only two highest degenerate bands (see Fig. [ 2 ]). The trial functions, which 
are used for the projection, were obtained from the digonalization of the site-diagonal den¬ 
sity matrix, calculated for the magnetically active bands in the basis of Ir 5 d orbitals.— > 2 ^ 
Namely, after the diagonalization of the density matrix, we pick up either 6 or 2 eigenstates 
(depending on the dimensionality of the model) with the largest eigenvalues and use them 
as the trial functions. Such construction guarantees that the Wannier functions are well 
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FIG. 3 . (Color online) Electronic structure of Sr2lr04 in LDA with the SO coupling, (a) Total 
and partial densities of states. The shaded area shows the contributions of the Ir Sd states. The 
positions of the main bands are indicated by symbols, (b) Band dispersion near the Fermi level, 
as obtained for the full LDA Hamiltonian in comparison the six-orbital models. The Fermi level is 
at zero energy (shown by dot-dashed line). 

localized in the real space: the main part of the density matrix with the largest eigenvalues 
is described by the “heads” of the Wannier functions, residing on the central site, and only 
small remaining part of this matrix is described by the “tails” of the Wannier functions, 
coming from the neighboring Ir sites. Thus, the main weight of the Wannier function is 
concentrated in its “head” part, while the contribution of “tails” is relatively small. Such 
procedure was extensively used in nonrelativistic calculations without the SO coupling.— 
The new aspect of the relativistic formulation is that the eigenstates of the density matrix 
become Kramers degenerate. Therefore, the trial functions and the Wannier functions {wi 
and W2) for each Kramer’s doublet can be chosen so to satisfy the conditions: \w2) = T\wi) 
and |tai) = —T\w 2 ), where T = i&yK is the time-reversal operation, written in terms of the 
spin Pauli matrix ay and the complex conjugation operator K. 

Then, the one-electron part of the model Hamiltonian ([T]) is identihed with the matrix 
elements of the LDA Hamiltonian in the Wannier basis: = {wia\i-LLDA\wjii). This 

procedure can be also reformulated as the downfolding of the LDA Hamiltonian.— >21 Then, 
the site-diagonal matrix elements describe the splitting of the atomic levels by the 
crystal field and the SO interaction, while the off-diagonal elements stand for interatomic 
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transfer integrals (or kinetic hoppings). 

The matrix of screened on-site interactions tj = \Uap^5\ has been calcnlated nsing sim- 
plihed version of the constrained random-phase approximation (RPA).— The RPA is nsed 
in the GW method in order to evalnate the momentnm and freqnency dependence of the 
screened Coulomb interaction, which is then used in the calculations of the self-energy.— The 
basic idea of constrained RPA is to switch off some contributions to the RPA polarization 
function (and, therefore, to the screening of U) related to the transition between the magnet¬ 
ically active bands (in our case, the Ir Sd bands).— The RPA is inadequate for this channel of 
screening (especially when it is evaluated starting from the PDA bandstructure) and should 
be replaced by a more rigorous method in the process of solution of the low-energy model 
(IT]) . The purpose of additional simplihcations is to replace the time-consuming RPA for the 
screening, caused by the relaxation of the atomic wavefunction and other (non-Sd) states, by 
much faster and more suitable for these purposes constrained LDA technique. After that, we 
consider (within RPA) the additional and most efficient channel of screening of the Coulomb 
interactions in the Ir 5 d bands by (the same) Ir 5 d states, which contribute to other parts 
of the electronic structure (mainly to the O 2 p and either Ba 5 d or Sr 4 d bands in Figs. [ 2 ] 
and | 3 |) due to the hybridization.— Such approximation incorporates the main channels of 
screening and, thus, reproduces reasonably well the values of static Coulomb interactions, 
obtained in full-scale constrained RPA calculations. The obtained matrix elements Ua/s^s 
have the following form: 

Uap^5 = j dr j dr'wl{r)wf}{r)vscv{TC,r')wl{r')ws{r'), (2) 

where the screened interaction nscr(i' 5 r 0 RPA is invariant under the time-reversal opera¬ 
tion and does not depend on the spin variables. 


IV. PSEUDOSPIN MODEL 

In this section we will consider the mapping of the electron model dl]) onto the magnetic 
model, formulated in terms of pseudospin variables Si = (iSf, iSf, iSf): 

7^5 = SjJ ijSj + SjgjH, ( 3 ) 

i>j i 
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where Jjj and Qi are the 3 x 3 tensors, describing interactions in the system of pseudospins 
and with the external magnetic held H, respectively. The pseudospin operators are repre¬ 


sented by the Pauli matrices: iSf = ^ 



qy _ I 


0 -i 
i 0 


, and S- = i 


1 0 
0 -1 


For each bond, Jij can be presented as the sum of its symmetric (S') and antisymmetric 
( 4 ) components: Jij = J\j + J\^ \ where = \{Jij + Jfj) and j\j = \{Jij — Jjj). 

^(S) 

The part Jb ^ incorporates all types of symmetric exchange interactions and its trace is the 
isotropic exchange interaction in the bond i-j: Jij = Tr Jb \ while Jb describes anisotropic 

^(A) 

DM interactions. Jb ^ has only three independent elements, which can be viewed as the 
components of some axial vectors (the so-called DM vector) = (db, db): 

/ 

"fM) _ 


0 

d^- 

% 

-dJ 

^ij 

\ 

—d^ 

0 

dx 


dJ- 

—d^- 

% 

0 

/ 


yielding the well known identity: SiJ^f'^Sj = dij[Si x Sj] 


A. Calculation of superexchange interactions 

In order to calculate the SE interactions, we adapt a standard procedure for the systems, 
whose degeneracy in the atomic limit is lifted by the crystal held and SO interaction. Namely, 
we assume that, in the atomic limit, the single hole resides on the highest Kramer’s doublet, 
obtained from the diagonalization of the site-diagonal part i = [tf£j] of the one-electron 
Hamiltonian. The states l^^i) and \(p2) of this Kramer’s doublet are used for the construction 
of eigenstates |±x), \±y), and |±2:) of the pseudospin operators and respectively. 

For convenience, we choose the phases of these states so that \ip2) = T\ipi) and \ipi) = 

-TM- 

Let us hrst explain the construction for |±^). For these purposes, one can choose any 
pair of states, which is obtained by the unitary transformation of \ipi) and |9?2)- Moreover, 
since the states are degenerate, the transformation will not change the total energy, and 
the model ([ 3 ]) will not contain the single-ion anisotropy term. Then, we employ the fact 
that, despite some complications caused by the strong SO coupling, the magnetic moment 
will always have a hnite spin component, and dehne the pseudospin states \+z) and \—z) as 
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those corresponding to, respectively, the maximal and minimal projections of spin onto the 
axis. The problem is equivalent to the diagonalization of the 2 x 2 spin Pauli matrix az in 
the basis of |(pi) and \^2)- 

Then, one can readily dehne two other groups of states as |±x) = -^\+z) ± and 

I±2/) = ±^|+z) + ^|—z). In this construction, the phases of |±z) were chosen to satisfy 
the condition: T\—z) = |+z) and T\+z) = —\—z). It allows us to define unambiguously all 
phases of \^z) but C, which transforms |±2;) to e^*^|±2;). The latter phase is dehned so to 
satisfy the condition {+x\ay\+x) = 0 . 


In order to hnd we evaluate the energy gain Tij{a, b), caused by the virtual excitations 
of the hole from the a-th orbital of the site i to the &-th orbital of the site j and vice versa, 
in the second order of perturbation theory with respect to the transfer integrals The 
denominators in the SE theory are given by the energies of charge excitations —>■ dfd^, 
which are the energies of the two-hole states. In the process of virtual excitations, the Pauli 
exclusion principle was guaranteed by the projection operators, which permit the hoppings 
only between occupied and unoccupied orbitals. Moreover, for the excited two-hole states, 
the problem was solved in the true many-body fashion, by finding the eigenstates and the 
eigenenergies from the diagonalization of the Coulomb interaction matrix U in the basis of 
= 15 Slater determinants, constructed from 6 atomic orbitals. This is a step beyond 
the mean-held approximation, which additionally stabilizes the AFM interactions.^^® Then, 
we consider all combinations of a and b = ±a:, ±y, or ±z, and map the obtained energy 
gains onto the pseudospin model ([ 3 ]) for iT = 0 . This procedure was discussed in details in 


Ref. 


30 


B. Calculation of g-tensor 


The ^f-tensor describes the interaction of the pseudopsin with the external magnetic held 


[see Eq. 


It can be evaluated using Eq. ( 31 . 34 ) of Ref. 


31 


from which one can hnd all 


9 elements of the tensor g at each site of the lattice: {+z\{Lx -|- ax)\+z) = {-\-z\{Ly + 

(^y)\+^) = {+z\{Lz + (Jz)\+z) = {-z\{Lx + ax)\+z) = g^^+ig^y, {-z\{Ly + ay)\+z) = 

gyx _j_ igvy^ and {—z\{Lz -|- crz)\+z) = g^^ ig^y, where L^, Ly, and Lz are the matrices of 
angular momentum operators in the Wannier basis. It is easy to separate the spin gs and 
orbital "g^ contributions to the ^f-tensor, by considering the matrix elements of only cr and 


11 




L, respectively. 


V. RESULTS AND DISCUSSIONS 

A. Two-orbital model for Ba 2 lr 04 

The two-orbital model is the simplest model, which can be considered. In Ba 2 lr 04 , the 
“j = 1/2” bands are separated from the rest of the spectrnm (see Fig. [2]) and the constrnction 
is rather straightforward. 

The form of transfer integrals in this case is very simple. Since t is hermitian, each 2x2 
matrix iij = [£“^] satishes the property: iji = ilj. Then, since all Ir sites are located in 
the inversion centers and connected by the translations, it holds tji = tij and, therefore, 
iij = ilj. Finally, since T-Llda is invariant nnder the time-reversal operation, we will have 
two more identities: (tjj)* = tfj and {tjj)* = —tfj, which can be obtained from = 

{Twia\T'HL]:)A\wjj^). Thus, in the two-orbital model, each Uj is proportional to the unity 
matrix t^j = tijl in the subspace spanned by the indices a{f5) = 1 and 2, where tij is a real 
constant. 

The behavior of tij is explained in Fig. 01 As expected, the strongest hopping occurs 
between nearest neighbors in the xy plane. There are also hnite hoppings between next- 
nearest neighbors in and between the planes. 

Since tij = tijl, all SE interactions in the two-orbital model are isotropic. They can be 
easily evaluated using the formula Jij = Atjj/U,^"^ where U = 1.52 eV is the effective on¬ 
site Coulomb repulsion, obtained in the constrained RPA for the two-orbital model. Then, 
using the values of transfer integrals, which are displayed in Fig. 01 we will obtain Jij = 
122.8, 2.5, and 0.8 meV for the nearest-neighbor (NN), next-NN, and interplane interactions, 
respectively. Since Jij > 0, all interactions are antiferromagnetic. 

B. Six-orbital model for Ba 2 lr 04 

The atomic t 2 g states are split into three doubly degenerate groups of levels, which in 
Ba 2 lr 04 are located at —209, —149, and 358 meV, relative to their center of gravity. Two 
lowest doublets correspond to j = 3/2, and the highest one - to j = 1/2. Thus, the 
splitting between the j = 1/2 and j = 3/2 states, which measures the strength of the SO 
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FIG. 4 . (Color online) Crystal structure and transfer integrals (in meV) associated with different 
Ir-Ir bonds in the two-orbital model for Ba2lr04. The Ir atoms are indicated by the big (red) 
spheres and the oxygen atoms are indicated by the small (yellow) spheres. 

coupling is very large. This justifies the use of the regular (nondegenerate) theory for the 
SE interactions. 

For the tetragonal compounds, the eigenstates \+z) (and \ —z) = —T\+z)), corresponding 
to the highest Kramer’s doublet, can be decomposed in the basis of xy, yz, zx, and x^—y"^ 
Wannier orbitals with both projection of spins: 

I + ^) = cly\Wxy,t) + + 4x\'Wzx,t) + cl2_y2\w^2_y2 + 

ciy\Wxy,i) + + Cl2_y2\w^^-y^i)- (4) 

Due to the symmetry constraint, the 3z‘^—r‘^ orbitals do not contribute to 1+^). The coef- 
hcients in this expansion depend on the relative strength of the crystal-held splitting and 
the SO interaction. They cannot be determined solely from the symmetry principles. For 
Ba 2 lr 04 , we obtain the following (nonvanishing) coefficients in the original lA/mmm coordi¬ 
nate frame: ci,„, = —fO.522, c!,^, = —ict, . = 0.603, and ci ,2 ,,2 = —0.004, which correspond 
to c^y = 0.004, cl^ = —ic^z = 0.426 -b iO.426, and c^ 2 _y 2 = —iO.522 in the lAi/acd frame. 
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FIG. 5. (Color online) Tensors of snperexchange interactions Jij (in meV), associated with different 


Ir-Ir bonds in the xy plane of Ba 2 lr 04 (top) and Sr 2 lr 04 (bottom). 


The Ir atoms are indicated 


by the big (red) spheres and the oxygen atoms are indicated by the small (yellow) spheres. 


For 


the sake of convenience, the parameters for both strnctnres are shown in the lAi/acd coordinate 


frame. 


The strongest transfer integrals, operating between the nearest neighbors in the xy plane, 
have the following form (in meV): 


/ -283 

0 

0 

±60 

0 

-776 \ 

0 

-283 

T60 

0 

-776 

0 

0 

T60 

-165 

0 

±792 

0 

±60 

0 

0 

-165 

0 

=F^92 

0 

i76 

=F*92 

0 

-226 

0 

^ *76 

0 

0 

±792 

0 

-226 ) 


where the npper (lower) sign stands for the bonds parallel to the x' {y') axis in the lA/mmm 
coordinate frame (see Fig. [S]). Here, the matrix is given in the local representation, which 
diagonalizes the site-diagonal part of the one-electron Hamiltonian [fi=j], as described in 
Sec. IIVAl Moreover, we adapt the following order of the Wannier orbitals: (m,s)= (1,1), 
(1, 2), (2,1), (2, 2), (3,1), and (3, 2), where m numbers the Kramer’s doublet in the increasing 
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FIG. 6 . The energies of two-hole states for Ba2lr04 (a) and Sr2lr04 (b), obtained using parameters 
of screened Coulomb interactions Uap-fS for the six-orbital model with and without the spin-orbit 
(SO) coupling. 

order of their energies and s number the states within each doublet. Similar to the 2- 
orbital model, the matrix elements of tij with same m do not depend on the s-indices and 
each such sub-block is proportional to the 2x2 unity matrix. However, there is a hnite 
coupling between states with different m’s. This coupling gives rises to the anisotropy of 
Jij. Moreover, since the signs of some of these matrix elements alternate between the bonds 
parallel to the x' and y' axes, the anisotropic part of Jij will also alternate in the x'y' 
plane. Another important factor, which is responsible for anisotropic properties of Jij is the 
intraatomic exchange interaction J It will be discussed below. Other parameters of the 
model Hamiltonian can be found elsewhere.— 

The form of the screened on-site interactions Uay-ys in the basis of relativistic Wan- 
nier orbitals is rather complex. Nevertheless, the main details of these interactions can 
be understood by considering the energies of two-hole excitations, which contribute to the 
SE processes (see Fig. |6]). These energies were calculated using the matrices of screened 
Coulomb interactions [Uay^ys], for which nscr(r,rO "'^ns obtained for two types of the elec¬ 
tronic structures: with and without the SO coupling [see Eq. ([2])]. In the case of perfect 
cubic environment and without the SO coupling, the two-hole states are split into three 
groups: ^Ti, degenerate ^T 2 and ^E, and ^Ai with the energies (W—SjT”), (U—J'), and 
{U+2J'), respectively,— in terms of the intraorbital Coulomb interaction U and the ex- 
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change interaction JT.— The tetragonal environment of the Ir^+ ions, realized in Ba 2 lr 04 , 
slightly lifts the degeneracy of the ^T 2 and states. The SO interaction further lifts the 
degeneracy of the ^Ti states. However, in all other respects the positions of the main energy 
levels in very similar with and without the SO interactions. The (averaged) parameters 
U and J can be evaluated from the centers of gravity of the three groups of levels. This 
yields: U = 2.86 (2.91) eV and J = 0.48 (0.49) eV with (without) SO interaction. Thus, U 
is generally larger in the 6-orbital model, in comparison with the 2-orbital one, due to the 
additionally screened by the j = 3/2 electrons, which is included in the 2-orbital model, but 
not in the 6-orbital one. 

The parameters of NN SE interactions in the xy plane are explained in Fig. [5l Since 
= J™ > these parameters favor the inplane configuration of the pseudospins, 
in agreement with the experiment.— Moreover, the phase of the off-diagonal element 
(in the Mi/acd coordinate frame) is bond-dependent, giving rise to the quantum compass 
interaction term. In the more conventional lA/mmm coordinate frame, the tensor Jij is 
diagonal with the parameters given by Jf-^' = Jf^ ± \Ji^\ and ip | where 

the upper (lower) sign stands for the bonds parallel to the x' {y') axis. The isotropic part 
Jij = = 109.8 meV is close to the value Jij = 123 meV, obtained in the 

2-orbital model. This is mainly because of the combination of two effects: On the one 
hand, U is larger in the 6-orbital model, which should lead to the smaller Jij. This decrease 
of Jij is partly compensated by somewhat stronger transfer integrals, operating between 
orbitals belonging to the highest Kramer’s doublet (—226 meV instead of —216 meV in the 
two-orbital model). 

As was already mensioned before, there are two important factors, which leads to the 
anisotropy of Jij\ (i) hnite transfer integrals, connecting the states with j = 3/2 and j = 1/2 
[see Eq. [3] and (ii) hnite intraatomic exchange coupling which lifts the main degeneracy 
of the virtual two-hole states (see Fig. [6]). For instance, using the same transfer integrals, 
but simplihed matrix of the screened on-site Coulomb interactions, which was reconstructed 
from the parameters of averaged U = 2.86 eV and JT” = 0, we have obtained totally isotropic 

o «■ 

tensor = Jjjl, where 1 is the 3x3 unity tensor and Jij = 71 meV. 

The direction of the uniaxial anisotropy is also controlled by the tetragonal crystal-held 
splitting /S.t 2 g between x'y' and doubly degenerate y'z' and z'x' orbitals without the SO 
coupling. If the x'y' orbital was located higher in energy, we would deal with the out-of- 
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plane configuration of pseudospins: Jf^ = . Such situation is indeed realized 

for Ajjg = 107 meV, associated entirely with the change of the hybridization due to the 
compression of the IrOg octahedra in the x'y' plane. However, the additional nonspherical 
Madelung interaction (see Ref. l25l) yields = —69 meV, thus changing the order of the 
t 2 g orbitals and enforcing the in-plane conhguration of the pseudospins > Jij), 

in agreement with the experimental data.— In mathematical terms, it leads to the inequality 


Cx'y'l < Wv'z'\ = Wz'x'\ ^he coefficients in Eq. 


y'z 


It is interesting that for the single-ion 


anisotropy (if the latter was appropriate in the analysis of some more general magnetic 
model), the preferential population of the y'z' and z'x' orbitals typically stabilizes the out- 
of-plane conhguration of spin and orbital magnetic moments. However, it should not be 
confused with the present situation, where we deal with the intersite interactions, which 
are governed by completely different processes rather than the single-ion anisotropy energy. 
This example emphasizes the importance of the tetragonal crystal-held splitting, which is 
sometimes ignored during the construction of the pseudospin models.— 

Due to the tetragonal M/mmm symmetry, the g-tensor of Ba 2 lr 04 has only two inequiv¬ 
alent matrix elements: and Other elements are identically equal to zero. The 

value of and g^^ are listed in Table H] together with the partial contributions of the spin 
and orbital components. 


TABLE I. Matrix elements of the g-tensor, obtained in the six-orbital model for Ba 2 lr 04 and 
Sr 2 lr 04 , and results of their decomposition into the spin (S) and orbital (L) parts (given in 
parenthesis). 


Compound 

^XX (ri^X ri^X\ 

y {ys '>yL ) 


9^^ {9l\97) 

Ba 2 lr 04 

1.796 (0.545,1.251) 

2.380 (0.909,1.470) 

0 (0,0) 

Sr 2 lr 04 

1.115 (0.208,0.907) 

3.332 (1.582,1.750) 

0.005 (0,0.005) 


C. Six-orbital model for Sr 2 lr 04 

In the case of Sr 2 lr 04 , the splitting of the t 2 g levels is —431, —4, and 435 meV. The 
symmetry properties of the |-l-z) orbital are given by the same Eq. (jl]) with the follow¬ 
ing (nonvanishing) coefficients: c^y = —0.015 =F f0.087, = —ic^^ = ±0.184 -|- iO.643, 
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and c^ 2 _y 2 = ±0.004 — zO.311, where the upper (lower) sign is referred to the site 1 (2), 
experiencing the counterclockwise (clockwise) rotation of the IrOe octahedra (see Fig. [T]). 

In the local representation, which diagonalizes the site-diagonal part of the one- 

electron Hamiltonian, the matrix of transfer integrals between sites 1 and 2 in the xy planes 
is given by (in meV) 


t 




( 218 + i60 0 0 t24tH 1 0 4 - i38 \ 

0 218-i60 ±24 tH 1 0 -4 - i38 0 

0 T24±ill -94-i69 0 ±29^*68 0 

±24±H1 0 0 -94 + i69 0 ±29 ± i68 

0 -4-i38 T29±i68 0 -144 - 0 

y 4 - i38 0 0 =f29 =f i68 0 -144 + J 


where the upper (lower) sign stands for the bond parallel to the x' {y') axis (see Fig. [T]for 
the notations). This matrix has both hermitian fh = 1(£^^. + tji) and antihermitian = 
— parts. The hermitian part has the same form as in Ba 2 lr 04 , where the off-diagonal 
matrix elements give rise to symmetric anisotropic interactions Jlj ■ The alternation of signs 
of some of these matrix elements will also lead to the alternation of anisotropic interactions 
in the xy plane. The antihermitian part is the new aspect, which is related to the fact that 
the neighboring Ir sites in the Mi/acd structure are no longer connected by the inversion 
operation. This part is responsible for the DM interactions. The transfer integrals, involving 
the highest Kramer’s doublet are generally smaller in Sr 2 lr 04 in comparison with Ba 2 lr 04 , 
mainly due to the additional rotation of the IrOg octahedra and deformation of the Ir-O-Ir 
bonds. Therefore, the SE interactions are also expected to be smaller in Sr 2 lr 04 . 

Due to the additional symmetry lowering, the matrix of the screened Coulomb interactions 
[Uay'ys] is even more complex than in Ba 2 lr 04 . Nevertheless, the energies of the two-hole 
states, obtained from [Uay^s], have the same “three-level” structure as in Ba 2 lr 04 , which 
is only slightly deformed by the lattice distortion and the SO interaction (see Fig. [6]). The 
averaged parameters U and JT” can be again evaluated from the splitting between these 
three groups of levels as U = 3.05 eV and J = 0.48 eV (both with and without the SO 
interaction). The value of J is comparable with the one in Ba 2 lr 04 . However, the Coulomb 
repulsion U is slightly larger in Sr 2 lr 04 . This behavior is consistent with the change of the 
electronic structure (see Figs. |2]and[3]): since the unoccupied Ba 5d states are closer to the 


18 




Fermi level and strongly hybridize with the Ir 5d states, the Conlomb U is expected to be 
more screened in Ba 2 lr 04 than in Sr 2 lr 04 .~ Moreover, it is reasonable to expect that the 
additional Mi/acd distortion in the case of Sr 2 lr 04 will make the t 2 g states more localized 
and, thns, the screening of U less efficient. This will fnrther reduce the values of the SE 
interactions in Sr 2 lr 04 . 

Considering only the values of interorbital Coulomb interactions W = U — ‘IJ = 1.90 
eV and 2.09 eV for Ba 2 lr 04 and Sr 2 lr 04 , respectively, we note a reasonable agreement with 


the results full-scale constrained RPA calculations reported in Ref. [m (U' is about 1.47 
eV and 1.77 eV for Ba 2 lr 04 and Sr 2 lr 04 , respectively). Moreover, the authors of Ref. [14 
used a simplihed M/mmm structure and theoretical lattice parameters both for Ba 2 lr 04 
and Sr 2 lr 04 , which may lead to the additional screening of U'. A more serious discrepancy 
is found for J7: our value of J7 is close to the atomic one, which seems to be reasonable, 
because J7 is only weakly screened in RPA.— However, the values of J7 reported in Ref. [14 
are about three times smaller, leading to the violation of the Kanamori rule W =IA — 2J7, 
presumably due to the contribution of the oxygen states to the Wannier functions.— This 
itself is an interesting point, because, according to Ref. Il3|, smaller value of J7 within the 
spherical model, which respects the Kanamori rule, should reduce the anisotropy of the 
exchange interactions. Therefore, it is interesting to which extent this anisotropy of the 
exchange interactions will be compensated by the anisotropy of the Coulomb interactions, 
which emerges in the full-scale constrained RPA calculations and manifested in the violation 
of the Kanamori rule. In any case, according to the analysis of the effective electron model 
based on the dynamical mean-held theory,— our values of the parameters U and W should 
correspond to the insulating behavior for Ba 2 lr 04 and Sr 2 lr 04 , thus justifying the use of 
the 1/U expansion for the analysis of the exchange interactions. 

The Mi/acd structure of Sr 2 lr 04 contains two Ir02 planes. The behavior of NN SE 
interactions in one of the plane is explained in Fig. |5l The parameters in another plane can 
be obtained by the 90°-rotation about the 2 ;-axis. As was expected, the isotropic part of the 
exchange interactions J 12 = +>^ 12 +'^ 12 ) ~ 39.0 meV is considerably smaller than in 

Ba 2 lr 04 . 

Since J 12 = Ju > J 121 ^^e pseudospins will favor the in-plane conhguration, similar to 
Ba 2 lr 04 and in agreement with the experimental sit nation. ->2^ In Sr 2 lr 04 , the parameter 
of the easy-plane anisotropy for the NN interactions, = 1 — J 12 /J 121 been recently 
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estimated in the X-ray resonant magnetic scattering experiments as 0.08,— which is close to 
our theoretical value of 0.087. The symmetric anisotropic part of J 12 is I = ^>^12 = 

0.73 meV, which is about one order of magnitude smaller than in Ba 2 lr 04 . This interaction 
is also bond-dependent. 

The antisymmetric part of J 12 can be represented in terms of the of DM vector (in meV): 
di 2 = (—0.1, —0.1, —3.97) (see Fig. [T]for the notations of the atomic sites). The phases of 
and alternate in the four NN bonds around the site 1. Therefore, since all NN atoms, 
surrounding the site 1, have the same direction of the pseudospin, the total contribution of 
and d^ to the canting of these pseudospins will vanish. On the other hand, the phases 
of d^ are the same for all NN bonds. Thus, d^ will be responsible for the ferromagnetic 
(FM) canting, which can be estimated as |df2/(2Jff)| ~ 2.8°. This value is smaller than 
the experimental estimate of 8°.- Nevertheless, the negative sign of d^ for the bond 1-2 is 
consistent with the counterclockwise rotation of the IrOg octahedra.— This picture can be 
also verihed experimentally.— 

The g-tensor relates the pseudospins with the value of true magnetic moments, which can 
be observed in the experiment. Using the value of = 1.115 (Table [I]), the local magnetic 
moment in the xy plane can be estimated as = 0.56 /ie, where the spin and orbital 
counterparts are ^gg^ = 0.10 /ib and \g^ = 0.46 pb, respectively. 


D. Calculations of Neel temperature 

Thus, the hrst-principles calculations have revealed a big difference of the magnetic mod¬ 
els in the case of Ba 2 lr 04 and Sr 2 lr 04 . On the one hand, the leading isotropic exchange 
interaction of 123 meV in Ba 2 lr 04 is about three times larger than that of 39 meV in 
Sr 2 lr 04 . In turn, the symmetric anisotropic interaction in Ba 2 lr 04 is an order of 

magnitude larger than in Sr 2 lr 04 . On the other hand, there is an appreciable DM interac¬ 
tion in Sr 2 lr 04 , but not in Ba 2 lr 04 . One of the puzzling points is that the experimental 
Neel temperature remains practically the same in both compounds (about 240 K). The aim 
of this section is to check whether such striking similarity can be explained using above 
parameters of interatomic exchange interactions derived in the SE approximation. 

Let us hrst investigate the effect of the DM interaction on the energy spectrum of the 
pseudospin model. In the Am/mmm coordinate frame, the exchange interaction tensor in 
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the bond 1-2, which is parallel to the y' axis, is given by 

/ 


J12 = 


- A Jr 


—J? 


0 


\ 




12 


V 


jzz 

^12 


>12 - u '12 

Jf2" + AJi2 0 
0 0 

where for simplicity we have dropped the small contributions of df 2 and d\ 2 - For the bonds 
parallel to the x' axis, A J 12 should be replaced by — A J 12 . By considering the transformation. 


J12 ^ J12 — UIJ12U2 


with 


u, = ul = 


\ 


COS( 

sine 

0 


sin0 0 
COS0 0 
0 1 




7 XX \ 
12 ) 


minimizes the energy of DM interactions,—!^ the tensor J 12 can 


•^12 = 


0 


V 


0 


JZZ 

'12 


and 0=1 arctan((if 2 / J] 
be transformed to 

^ 7 ^ 2 " - A J 12 

0 Jf2" + AJi2 0 
0 0 J 

where Jf|^ = Jf 2 *\/l + (<^i 2 /'^f 2 )^- Thus, the DM interactions alone do not conhne the 
pseudospins in any particular directions.—!^ Moreover, after such transformation to the local 
coordinate frame, the effect of the DM interactions can be combined with Jf^- Since in the 
6-orbital model for Sr 2 lr 04 , d \2 = 3.97 meV while J^f = 40.2 meV, the renormalization of 
J 12 due to the DM interaction is only about 0.5 %. Therefore, we conclude that the effect of 
the DM interaction on the energy spectrum is small and can be neglected and, as far as the 
energy spectrum is concerned, the main ingredients of the pseudospin model are essentially 
the same in the case of Ba 2 lr 04 and Sr 2 lr 04 . 

Below, we will concentrate on two mechanisms of the magnetic ordering in iridates: the 
hrst one is due to the in-plane anisotropy, which emerges in the 6-orbital model, and the 
second one is due to the interlayer exchange coupling, which is relevant to the 2-orbital 
model of Ba 2 lr 04 . Thus, we consider the following general compass Heisenberg model 

E ■SfSI + i 5 : (J||5f5J + Jx5?sj) 


<ij> in plane <p>||a: 


1 E + J 




( 6 ) 


<ii>\\y 


<ij> inter plane 
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where it is convenient to introduce the shorthand notations: Jz = Jf|, J\\ = Jff + AJ 12 , 
J_L = Jf 2 — AJ 12 , and J' is the coupling between the atoms, which belong to different 
planes, separated by the primitive translation c along the z axis. The magnon spectrum of 


this model for J' = 0 was calculated in Ref. 


42. It reads 


(4Jav + Bq — 74q)(4Jav + Bq + Aq + Javi?), 

(4Jav -B,-A, + Jav^7)(4Jav -B, + A,), 

where Jav = {J\\ + J±)/2 = 


( 7 ) 


= {J± + Jz) COSqz: + (J|| + Jz) COSQy, 

Bq = {J± - Jz) cosqz, + (J|| - Jz) cos qy, (8) 


and 

g = 0.16(J|| - J±)y{JlS) = 0M{AJ,2/J^^)VS (9) 

is the quantum gap. Moreover, we have introduced the renormalization factor ^ = 1 + 
0.0785/S', which is taken equal to its value in the two-dimensional Heisenberg model. Then, 
for small q we obtain 

4Jav ~ Bq — Aq —>■ J±q‘^ “t" 

4Tav — Bq + Aq —)■ 4(Jav “t" J z), 

4^av + Bq — Aq —)■ 4(Jav — Jz) + Jz(f = Jz ijf + /) ) 

4^av Bq + Aq ^ 8Jav, (10) 


where the parameter / describing the in-plane symmetric anisotropy is dehned as 


/ = 4(Jav - 


( 11 ) 


Therefore, we have: 

^ ^CV8JavJ.(g2 + /), 

^ SC^J^{J.. + Jz){J±ql + J\\ql + J..g). ( 12 ) 

The hrst mode is related to the out-of-plane pseudospin rotation, while the second corre¬ 
sponds to the in-plane rotation. 

To obtain magnetic transition temperatures, we map the Heisenberg model ([6]) onto 
the non-linear sigma model, having the same excitation spectrum, Eq. oa. see Appendix. 
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Treating the magnetic excitations, as slightly different from the case of the XY anisotropy,— 
we obtain in the regime / S> ma.x{a,g), a = 2J'/J, the following eqnation for the Neel 
temperature (see Appendix): 

+ __ 

^op^ip fr In^ (// m&x{a,g)) 

where A ~ 3.5, Cop = y/SJ^.vJzSC and Qp = ■Y/4Jav(Jav + Jz)SC are the out-of-plane and 
in-plane spin-wave velocities, Ps = + l/pav)~^ is the effective spin stiffness (pz,av = 

Jz,a.vCSSo), fr = f{So/S)^ is the renormalized anisotropy parameter, Sq = 0.303 for S = 1/2 
is the ground-state magnetization. In the absence of compass anisotropy, / = 5^ = 0, we 
obtain instead^ 

Tn = 4:71 ps I In - hdln^^-o.oel , (14) 

CopCjpO^T- -/tv J 

where = q:(S'o/S') is the renormalized interlayer coupling parameter. 

The parameters and the resulting magnetic transition temperatures are listed in Table [TTl 
Lets us first discuss the results of the 6-orbital models for the Ba2lr04 and Sr2lr04. Judging 

TABLE II. Parameters used in Eqs. (| 13 p and (| 14 p for the transition temperature equation and the 
calculated T/v in different regimes (the values of Jav, J\\ — J^_■, and Jz are in meV, T/v is in Kelvins, 
and other parameters are dimensionless). 



'^av '^11 

- J± 

Jz 

/ 

9 

a 

^ ^a=0,g=0 

^op ^ip J- TV 

Tn 

Ba2lr04 (2-orb.) 122.8 

0 

122.8 

0 

0 

1.4-10-^ 

200.9 200.9 

239 

Ba2lr04 (6-orb.) 

110.6 

15 

108.1 0.09 6-10-3 

1.3-10-^ 

178.9 179.9 371 

414 

Sr2lr04 (6-orb.) 

40.2 

1.5 

36.7 

0.38 4-10-^ 

1.5-10-3 

62.8 64.3 181 

216 


Tn = ^TTPs < In 



by the ratio between the anisotropy parameters /, g and interlayer isotropic parameter, a, 
we have the relation f ^ g ^ a, which holds for both compounds. Thus, the in-plane 
anisotropy is expected to be mainly responsible for the magnetic ordering. The differences 
between in-plane and out-of-plane components of the symmetric anisotropy tensor, (Jav —<Jz), 
are close to each other and equal to 2.5 meV (in Ba2lr04) and 3.5 meV (in Sr2lr04). 
However, due to the difference in the absolute value of J^, we obtain completely different 
anisotropy parameters /, g and, therefore, the transition temperatures. For Sr2lr04, the 
calculated temperature of 216 K is in the good agreement with the experimental value of 
240 K. This is consistent with the finding of Jackeli and Khaliullin (Ref. Il3l) . who used 


23 















the experimental Tat in order to estimate the values of the exchange interactions and these 
values are close to ours. However, the situation is different in the case of Ba 2 lr 04 , where 
the theoretical T/v is overestimates by factor two. Interestingly, in the case of the 2-orbital 
model for Ba 2 lr 04 , which, in analogy with the cuprates,— contains only in-plane and inter¬ 
plane isotropic exchange interactions, we observe a good agreement between theory and 
experiment. However, this agreement is probably fortuitous. 


VI. BEYOND SUPEREXCHANGE 


The main purpose of this section is to discuss the effect, which are not included to the 
regular SE model. Our main concern is the following: since the SE model is based on 
the second-order perturbation theory for the transfer integrals, it implies that all effects 
of the SO coupling, which are included to these transfer integrals, are also automatically 
treated only up to the second order. Since the SO coupling is large in iridates, this may be 
rather crude approximation, which does not take into account some important anisotropic 
interactions. For instance, in the mean-held approximation for the SE model, all pseudospins 
in the single xy plane of Ba 2 lr 04 and Sr 2 lr 04 can rotate rigidly at no energy cost. Besides 
quantum effects, considered in the previous section (see also Ref. Il6|), this maybe related to 
the lack of the in-plane anisotropy, which typically appears only in the fourth order of the 
SO coupling. 

If we wanted to include these effects in the pseudospin model ([3]), our strategy would be 
to go beyond the second order perturbation theory for the transfer integrals and consider 
higher-order terms, which give rise to the new type of interactions, such as biquadratic and 
ring exchange.—*^ They will affect both anisotropic and isotropic parts of the exchange 
interactions. Therefore, in the pseudospin formulation, based on the strong SO coupling, 
these two types of the effects are connected with each other: if we want to consider the 
higher order anisotropic interactions, we have to deal with the biquadratic and ring exchange 
terms, which will affect all other exchange interactions, including the isotropic ones. Such 
pseudopsin Hamiltonian is no longer presented in the bilinear form (j3]). 

Nevertheless, in the present work we take a different strategy and in order to evaluate the 
higher-order contributions (and, therefore, check the validity of the SE model) in Ba 2 lr 04 
and Sr 2 lr 04 , we solve the original electron model ([1]) in the mean-held HE approximation. 
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FIG. 7. (Color online) Results of constrained unrestricted Hartree-Fock calculations for Ba2lr04 
in the staggered “antiferromagnetic” field = 0.68 meV. The direction of the magnetic held in 
the planes z = 0 and z = 1/2 is specihed by azimuthal angles {(p+Ap) and {(p—Ap), respectively 
(in the lAi/acd coordinate frame), (a) is the total energy dependence on cp for Acp = 0. (b) is 
the total energy dependence on Acp for cp = 0. (c) and (d) explains the geometry of the staggered 
magnetic held for (a) and (b), respectively. 


where we also apply the staggered external magnetic field, which controls the directions of 
the spin and orbital moments. We have fonnd that the held of (IbH = 0.68 meV is generally 
sufficient for these purposes. 

The weak point of the HF approach is that it treats all on-site electron-electron interac¬ 
tions on the mean-held level, whereas in the SE theory such processes are treated rigorously 
by solving the exact eigenstates problem for the virtual two-hole states. However, in this 
particular case, we do not expect large error caused by the mean-held approximation (some 
comparison for transition-metal perovskite oxides can be found in Ref. l25l) . On the other 
hand, the HF method does not employ any additional approximations regarding the relative 
strength of transfer integrals and the on-site Coulomb repulsion and, in this sense, is the 
more superior approach in comparison with the SE theory. 

Let us start with Ba2lr04. The geometry of the constraining held in this case is explained 
in Fig. [71 First, let us consider the case, where the helds in the two adjacent planes z = 0 
and z = 1/2 are rotated in phase. Then, the total energy exhibits the minimum at 0 = 0 
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FIG. 8. (Color online) Tensors of superexchange interactions Jij (in meV and in the I^\/acd 
coordinate frame), associated with different Ir-Ir bonds between adjacent planes in Ba 2 lr 04 . The 
Ir atoms are indicated by the big (red) spheres and the oxygen atoms are indicated by the small 
(yellow) spheres. 


(modulo TT, in the Mi/acd coordinate frame). This effect can be actually included in the 
SE model and is related to the anisotropy of the exchange interactions between adjacent 
planes.— The behavior of these interactions is explained in Fig. [HI Then, the mean-held 
energy of the magnetic order, depicted in Fig. [3::, is given by F(0) = — A Jout cos 20 (per 
one Ir site), where AJ^ut = l^out “ '^outl- 

Thus, in the SE approximation, the energy should remain invariant with respect to the 
antiphase rotations of the pseudospin (Fig. [711). In the other words, if we hx 0 and consider 
the conhgurations, where the directions of the pseudospins in the adjacent planes z = 0 
and 2 ; = 1/2 are specihed by the azimuthal angles (01-A0) and (0—A0), respectively, 
the mean-held energy of such conhgurations should not depend on A0.— This property is 
indeed strictly observed when we use the exchange parameters, derived in the SE model. 
Because of this degeneracy, the authors of Ref. [l6j had to go beyond the mean-held theory 
and consider the ehect of the quantum huctuations in order to explain the experimentally 
observed magnetic ground-state structure of Ba 2 lr 04 (corresponding to 0 = A0 = 0 in 
the lAi/acd coordinate frame).— The most interesting aspect of our analysis is that this 
degeneracy can be lifted even on the mean-held level if one goes beyond the SE model 
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and consider more rigoronsly the higher-order contribntions of the transfer integrals in the 
framework of the unr estricted HF calcnlations. The dependence of the HF total energy 
on A0 is shown in Fig. [H] (for 0 = 0). It clearly shows that the higher order-anisotropic 
interactions, which are inclnded in the HF calcnlations, lifts the degeneracy and stabilizes 
the experimentally observed magnetic gronnd state. The energy barrier, cansed by these 
interactions, is abont 4.5 peV, which is at least comparable with the effect of qnantnm 


f. 3. 


flnctnations considered in Ref. Il6|. Thns, the effect is robnst and cannot be neglected in the 
realistic analysis of the magnetic properties of Ba 2 lr 04 . 

Next, we evalnate the effect of biqnadratic exchange on the NN interaction Jf| in the 
a:|/-plane of Ba 2 lr 04 . If the magnetic properties of some material were indeed described 
by the bilinear Hamiltonian ([3]), the valnes of the exchange parameters wonld not depend 
on the method, which is nsed for their calcnlations. For instance, in the mean-held HF 
method, one could evaluate Jf| from the total energy difference between FM and AFM 
states, by aligning the magnetic moments parallel to the z axis: Then, 

if the bilinear parametrization ([3]) for the magnetic Hamiltonian were indeed appropriate, 
this value of Jf| should be close to the one obtained in the SE model. Nevertheless, the 
straightforward HF calculations yield Ey^—E^^ = 83.8 meV, which is 22% smaller than 
Jfl = 108.1 meV, obtained in the SE model. This deviation is the measure of biquadratic 
(or ring-type) exchange interactions, existing in the system. Thus, as expected from the 
discussion in the beginning of this Section, the higher order anisotropic effects in Ba 2 lr 04 
coexist with appreciable biquadratic contributions to the isotropic exchange interaction. In 
this sense, we obtain very consistent description for Ba 2 lr 04 . Unfortunately, we could not 
obtain a stable in-plane FM solution in the HF method and, thus, evaluate the in-plane 
elements of the exchange tensor from the total energy difference. Generally, one can expect 
similar contribution of biquadratic interactions to the in-plane and out-of-plane components 
of the exchange tensor. 

The behavior of Sr 2 lr 04 appears to be rather different from Ba 2 lr 04 . Since the transfer 
integrals are smaller in Sr 2 lr 04 , while the Coulomb interactions are slightly larger, it is 
reasonable to expect that the t 2 g states are more localized in Sr 2 lr 04 , which additionally 
justihes the use of the SE model. This is indeed what we have obtained by comparing results 
of HF calculations and the SE model. The fact that Ba 2 lr 04 appears to be “more itinerant” 
than Sr 2 lr 04 can be seen already from the comparison of the band gap, obtained in the 
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HF method for the AFM ground state, which is substantially smaller in Ba 2 lr 04 (1.3 eV, 
against 1.8 eV in Sr 2 lr 04 ). It should be noted, however, that the HF gap is considerably 
larger than the experimental one, due to the lack of quantum and thermal fluctuations, as 
was conhrmed by the DMFT calculations.— 

First, we consider the HF solutions for the FM and AFM states, where all magnetic 
moments are parallel to the ^ axis. The total energy difference between these states is 
31.7 meV, which is much closer to the value Jf| = 36.7 meV, obtained in the SE model 
(the difference is about 14%, which can be again regarded as the measure of biquadratic 
interactions in the system). In Sr 2 lr 04 , it is practically impossible to evaluate the in¬ 
plane elements of the exchange tensor from the total energy difference: because of the DM 
interaction, the in-plane FM state is unstable and converges to the AFM state (with small 
FM canting of the magnetic moments). 

Next, we consider the higher-order anisotropy effects in Sr 2 lr 04 . For these purposes we 
take the weakly FM state and rotate magnetic moments by the external magnetic held of 
/ieFT = 0.68 meV, which couples to the weak FM moment in the xy plane. The results of 
such constrained HF calculations are summarized in Fig. |9l We note the following: (i) The 
total energy depends on the direction of the magnetic moments in the xy plane. However, 
this dependence is very weak (the characteristic energy barrier is about 0.25 meV, which 
is an order of magnitude smaller than in Ba 2 lr 04 ); (ii) The angle (A0) between magnetic 
moments of the sites 2 and 1 (see Fig. [1] for the notations) is nearly constant, meaning 
that it is mainly controlled by the DM interaction dl 2 , while the effect of other anisotropic 
interactions (that are not taken into account in the SE model) are relatively small. Since the 
energy gain caused by the DM interaction is proportional to df 2 sin A(j), the obtained values 
of —270° < A0 < —180° are well consistent with the sign d \2 < 0 of DM interactions for 
the counterclockwise rotation of the IrOg octahedra around the site 1 (see Fig. [1]). Yet, one 
interesting aspect of the HF analysis is that the angle A0 is different between, separately, 
spin and orbital magnetic moments. Without external held {H = 0), A0 is about —185.2°. 
It corresponds to the FM canting of 2.6°, which is close to 2.8°, obtained in the SE model. 
The values of spin and orbital magnetic moments, obtained for the in-plane (out-of-plane) 
magnetic alignment are 0.13 and 0.48 /ib (0.71 and 0.83 /xb), respectively, which are in good 
agreement with the values of corresponding matrix elements of the g-tensor, reported in 
Table H] for the SE model. 



FIG. 9. Results of constrained unrestricted Hartree-Fock calculations for Sr 2 lr 04 in the magnetic 
field = 0.68 meV, which couples to the weak ferromagnetic moment in the xy plane. The 

direction of the field is specified by the azimuthal angle (p. The upper panel displays the behavior 
of the total energy: the symbols show calculated points, while the solid line is the result of interpo¬ 
lation E{4)) = ABcos4:(p. The lower panel shows the angle between spin (Aps)^ orbital (Api) 
and total {Ap) magnetic moments of the sites 2 and 1 in Fig. [TJ 

Thus, we obtain a very consistent description also for Sr 2 lr 04 : (i) To a good approxima¬ 
tion, the magnetic Hamiltonian has the bilinear form (I3l), inherent to the SE model; (ii) The 
higher-order anisotropy effects, beyond the SE model, are negligibly small. This makes the 
main difference from Ba 2 lr 04 , where (i) the deviations from the bilinear form are significant 
and (ii) the higher-order anisotropic exchange interactions are important. 


VII. SUMMARY AND CONCLUSIONS 

The main purpose of this work was to critically evaluate the abilities of the SE model 
for the analysis of magnetic properties of the layered iridates Ba 2 lr 04 and Sr 2 lr 04 . Being 
based on the hrst-principles electronic structure calculations with the SO coupling, we have 
first derived the effective low-energy electron model for the t 2 g bands, which are located 
near the Fermi level and primarily responsible for the magnetic properties of Ba 2 lr 04 and 
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Sr 2 lr 04 . This electron model was further mapped on the pseudospin model using the theory 
of SE interactions in the limit of large on-site Coulomb repulsion. We have clarihed the 
microscopic origin of the bond-dependent anisotropic exchange interactions, as well as the 
antisymmetric DM interactions, caused by the anti-phase rotations of the IrOe octahedra in 
Sr 2 lr 04 . The pseudospin Hamiltonian problem has been solved by means of the non-linear 
sigma model, that has hnally allowed to evaluate the Neel temperature for both considered 
compounds. We have demonstrated that while for Sr 2 lr 04 the theoretical Neel temperature 
is in good agreement with the experimental data, for Ba 2 lr 04 it is overestimated by factor 
two. We have argued that this discrepancy is quite consistent with the limitations of the SE 
model for Ba 2 lr 04 , which is the more “itinerant” system than Sr 2 lr 04 . Such “itineracy” is 
directly related to the details of the electronic structure of Ba 2 lr 04 : the lack of rotations of 
the IrOe octahedra and the proximity of the Ba 5d states to the Eermi level make the t 2 g 
bandwidth increase and more efficiently screen the Coulomb interactions in this band. Thus, 
the t/U expansion for the magnetic energy converges slower and higher-order terms, beyond 
the SE contributions, start to play an important role. Since the effect of SO interaction in 
the SE formulation is included to the transfer integrals, the higher-order terms automatically 
improve the description also for the anisotropic exchange interactions. In fact, by solving the 
low-energy electron model for Ba 2 lr 04 in the HE approximation, we were able to reproduce 
the experimental magnetic ground states structure of this compound even on the mean-held 
level, without invoking to quantum effects. 
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Appendix: Derivation of the non-linear sigma model for compass Heisenberg model 
and its renormalization 

1. Nonlinear-sigma model 

To obtain the action of the continuum model we pass to the coherent state representation 
for spin operators and represent the corresponding vectors of spin directions following the 
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standard procedure 


5, = (-ir^n,Vl - (L./^)2 + L,, (A.l) 

where Lj ■ iij = 0, = 1, and the helds n* and Lj represent the staggered and uniform 

component. Substituting Eq. flA.lIl into (E]) we obtain the Lagrangian: 

T c2 _ c2 

^[n, L] = —“ Y ('^ll«+<5. + 

'll X 

+ 2 (jzLm+s + {.Jzn\ + ■h'^^s + —Y^} 

1^5 ^ ' 

+\'£(-hii’tLu.+ 

+Y Y “ s LrYxa^iii], (A.2) 

i,6y i 

where = nx,y,n-^ = Uy^x, and we keep only terms, which do not vanish and give the 
leading contribution in the continuum limit. Expanding 

ni +5 = n* + (5V)n. + ^{S'^S^dadb)n. + ... (A.3) 

and similarly for Lj+ 5 , we obtain: 

5*2 r 

L[n, L] = Y / dY [J^{Vn,f + J\\{dxnxf + JwidynyY + J±{dxny f + J±{dynxf + Jzfnl] 
+ - y cfx. [2{AJz + (J|| + J± — 2Jz){nl + ny))L‘f+2{J\\ + J±){Ll + Ly)] 

+ iL ■ [n X (^-4) 


where we have dehned / according to flTT|) . Performing the integration over L, we hnd 


5*2 r 

L[n] = Y [Jz{Vn^f + J\\{dxnxf + JuidyHy^ + J±{dxnyy + J±{dynxf 


1 

2 

1 

2 


Yx 

Yx 


^ f d"x 


2(4J, + (J|| + Y - 2J,)(Y + Y)) 

2(4J, + (J|| + Y - 2J,)(1 + Y + Y)) 

1 r ^ 12 

2(4J, + (J|| + - 2J,)(1 +nl + nl)) 


(A.6) 
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In the following we assnme the preferable direction of magnetic order along the y axis. 
Representing 


and expanding in nx,z we obtain two branches of the magnon spectrum 


El = 4SV.(J|| + Jx) (?" + /) , 

El = 2S'‘{2J. + J, + Ji)(J||g2 + 


which coincides with small q expansion of the results of Sec. IV PI and Ref. 
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(A.6) 


(A.7) 


2. Perturbation theory 


In the following we concentrate on the classical part of the Lagrangian (IA.5I1 . renormalized 
by the quantum fluctuations, 


Lci[n] = - y {pr [{Vn^f + + p\\ [{d^rixf + [dyUyf] + p± [{d^Uyf + [dytixf] } . 

(A.8) 


In Eq. flA.Sp we use the quantum-renormalized spin stiffnesses, pr = JzSSqC, and p\\^±_ = 
J|l_j_S'S'oC, where S'o = S' — 0.197 is the ground state magnetization of the square-lattice 
Heisenberg model, C, = 1-1-0.157/(25') is the exchange parameter renormalization factor, the 
bare spin stiffnesses anisotropy, and the renormalized easy plane anisotropy fr = /5q/(S'/')^. 
Following the standard procedure^*^, we assume that the excitations, described by Lci[n] 
are cut on the ultriviolet at the momentum Auv = R/c, where c = \/8JSC is the renor¬ 
malized spin-wave velocity; the remaining (non-universal) contribution of the other part of 
momentum space yields the abovementioned quantum renormalizations. 

Assuming again the long-range order along the ?/-axis, introducing tt = [rix, Uz), and using 
Eq. flA.bp . we obtain 

1 


LcJnl = 


Pr [(Vvr^)^ frT^‘f\ + P\\{dxTlxf + P±{dyTlxf 


(A.9) 


vr^ 


TT^ 


■^^^{’Kdy'Kf + \ + ^ I '^^xln(l - - h / d^x\/r- ^2 


vr^ 


where the hrst term in the second line comes from the integration measure and in the last 
term we have introduced external magnetic held along y axis, which will be put to zero in 
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the end. To perform renormalization of Eq. flA.9D . we decouple the interactions via the 
Wick theorem 


Lci[n] = - y [pr [(VtT;,)^ + /^TT^] + + p^idyir^f + p\\{Til){dyT:af 

+pdT^l){d^'naf + p\\{{dy'Kaf)TTl + p^{{d^'KafWa\ 

^ ^ d^x [hTT^ + ■^(3(7r2) + + ^(3(7r2) + -Ttt^ 


(A.IO) 


Rescaling the helds 


t ZxT^xi 
TT^ ^ Z^T^zi 


we obtain renormalized parameters: 


Pr — \j)x + Pxyizz)\ 1 
P\\R = Zl [pii +P±(vr^)] , 
P±R = Zl [p± + P||(vr^)] , 


PrJr + ha — Z‘\ h -\- prfr + T 

(fq 


(Pq 


P\\<fy+ P±<ll 


(27r)2pr(g^ + fr) + h 


Pxy,RgR + hn — Zx 


h + T 


P\\(fy+ P±(ll 


(2vr)2 p\\ql + Pxg 2 ^ ^ 


T ++ {t^D) , 

-^+^(3(vr^) + , 

(A.ll) 


where pxy = (py +p±)/2, 


= T 
(tTx^) = T 


(Pq 1 

{ 27 iP Pr{q^ + /r) + h’ 

jPq _1_ 

(27r)2 piiql + p^ql + PxyQ + h’ 


(A.12) 


and qr is the gap, generated for tTx mode, which also contains the non-universal bare value 
p, determined by the equation ([9]). From the equations flA.lip we obtain 


Zx — Zz — Z, 

Hr = z‘^h 1 + + (^x)) 

PrIr = PrZ^fr [l - (vr^)] , 

Pxy,RgR — Pxy,rZ P [l ("^x)] ' 


(A.13) 
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Finally, Z is fixed by the condition, that renormalizes the same way, as Tr^,, which is due 
to 90° rotation symmetry in the plane. This implies Hr = Zh, such that 


Z = 1 - i ((x^) + (xS), 

PR ^ Pr[l- (ttI)] , 

Jr ~ /r [1 - 2(7r2)] , 

Pxy,R/Pxy = 1 ~ z) t 

7 i ?/7 = = 1 - 2(vr^), (A.14) 

where we have introduced 7 = (py — Pa_)/{‘ 2 ^pxy) and neglected small anisotropy terms in the 
second and third lines. Being rewritten through these quantities, the effective Lagrangian 
reads 

iRl"] j + frp] + Pxy.R [(Vlj)^ + flS + 

+ YZ? ] ■ (A-15) 


3. Renormalization group 


To perform RG analysis we introduce sharp momentum cutoff at scale A and vary A from 
Auv to the smallest possible scale; in the following we replace accordingly the index R at the 
renormalized quantities by A, denoting explicitly, at which infrared scale they are evaluated. 
We also assume in the following that f > g > a. According to the obtained expressions, we 
perform renormalization group procedure in two steps. At the hrst step we integrate degrees 
of freedom at momenta scales // < A < A^v In this range we can neglect small difference 
between x— and 2 ;— modes in equations flA.141) and obtain the standard flow equations of 
the 0(3) non-linear sigma model with small easy-plane anisotropy 


dt\ 

(iln(l/A) 
dln^A 
(iln(l/A) 
din Pa 
dln(l/A) 


J.2 I J.3 


— tR, 


dlnyA 

dln(l/A) 


din /a _ 
dln(l/A) 


(A.16) 
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where = T/{27ip\). In the first Eq. of flA.lhp we have added the two-loop term of the 
0(3) model. The solution of Eqs. flA.lhp is well known, 







9a 

9r 7 


1 — In 



4 = (hiX 

fr \tAj 


1 



(A.17) 


where we have introduced U = T/{27ipr). The scaling is stoped at A = A/, which fulhlls 
/aj. = Aj. The condition tA/ = 1 determines the Kosterlitz-Thouless temperature in the 
absence of the in-plane anisotropy. At the scales < A < the ^-mode is fully gaped, 
and we obtain behavior of the coupling constants 

dtxy,A _ din /a _ 

(iln(l/A) dln(l/A) ’ 


dln^A 

(iln(l/A) 

dln^A 

(iln(l/A) 


"Oy,A/2, 


din 7 a _ 

(iln(l/A) 


(A.18) 


which is in the XY universality class. The consideration of this regime is similar to the case 
of quasi-two-dimensional easy plane model,— and yields the result for the Neel temperature 
in Eq. (IT^ of the main text. 
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